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INTRODUCTION 

An increasing number of studies advocate the view that cere- 
bellar contributions to behavior are not confined to motor 
control but also extend to cognitive function (e.g., Stoodley 
and Schmahmann, 2010). Consistent with this, convergent 
evidence from clinical, neuroimaging and anatomical trac- 
ing studies in primates suggests that the cerebellum forms 
"closed-loop" connections with neocortical brain regions includ- 
ing the prefrontal cortex (Middleton and Strick, 1994, 2001; 
Kelly and Strick, 2003; Schmahmann, 2004; Allen et al., 
2005; Krienen and Buckner, 2009; Strick et al, 2009; O'Reilly 
et al., 2010; Stoodley and Schmahmann, 2010; Buckner 
et al., 2011; Stoodley, 2012). These anatomical connec- 
tions provide the neural basis through which cerebellar 
contributions to neocortical processing may occur, enabling 
integration of sensorimotor information across hind- and 
fore-brain. 

The understanding of such distributed networks may be espe- 
cially pertinent given that abnormal prefrontal-cerebellar interac- 
tions are implicated in disorders such as autism and schizophrenia 
(Andreasen et al., 1996; Andreasen and Pierson, 2008; Fatemi 
et al., 2012). In particular, imaging studies frequently report 
abnormalities of the cerebellar vermis in schizophrenia (e.g., 



Although recent neuroanatomical evidence has demonstrated closed-loop connectivity 
between prefrontal cortex and the cerebellum, the physiology of cerebello-cerebral circuits 
and the extent to which cerebellar output modulates neuronal activity in neocortex 
during behavior remain relatively unexplored. We show that electrical stimulation of the 
contralateral cerebellar fastigial nucleus (FN) in awake, behaving rats evokes distinct 
local field potential (LFP) responses (onset latency ~13ms) in the prelimbic (PrL) 
subdivision of the medial prefrontal cortex. Trains of FN stimulation evoke heterogeneous 
patterns of response in putative pyramidal cells in frontal and prefrontal regions in both 
urethane-anesthetized and awake, behaving rats. However, the majority of cells showed 
decreased firing rates during stimulation and subsequent rebound increases; more than 
90% of cells showed significant changes in response. Simultaneous recording of on-going 
LFP activity from FN and PrL while rats were at rest or actively exploring an open 
field arena revealed significant network coherence restricted to the theta frequency 
range (5-10 Hz). Granger causality analysis indicated that this coherence was significantly 
directed from cerebellum to PrL during active locomotion. Our results demonstrate the 
presence of a cerebello-prefrontal pathway in rat and reveal behaviorally dependent 
coordinated network activity between the two structures, which could facilitate transfer 
of sensorimotor information into ongoing neocortical processing during goal directed 
behaviors. 
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Okugawa et al., 2007; Lawyer et al., 2009; Henze et al., 2011) 
and direct electrical and transcranial magnetic stimulation of the 
vermis has shown some efficacy in treating the cognitive and emo- 
tional symptoms of the disease (Heath, 1977; Demirtas-Tatlidede 
etal., 2010). 

Although the presence of a cerebello-cortical reciprocal net- 
work has not been demonstrated in non-primates, electrophysio- 
logical and amperometric studies have highlighted, respectively, 
the existence of a prefrontal-olivo-cerebellar pathway in anes- 
thetized rats (specifically to vermal lobule VII; Watson et al., 
2009), and modulation of prefrontal dopamine release follow- 
ing cerebellar stimulation in anesthetized mice (Mittleman et al., 
2008). Anatomical data also suggest the existence of disynap- 
tic fronto-cerebellar connectivity in rat (Suzuki et al., 2012) and 
preliminary data obtained in mouse suggest a neural connec- 
tion exists between the cerebellar nuclei and prefrontal cortex 
(Arguello et al., 2012). Recent evidence has also highlighted 
the importance of cerebellar plasticity in goal-directed behav- 
ior and spatial navigation in mice (Burguiere et al., 2010; 
Rochefort et al., 2011). Together, these studies suggest the basis 
of a rodent prefrontal-cerebellar network reminiscent of that 
described anatomically in primates (Middleton and Strick, 2001; 
Kelly and Strick, 2003). 
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While the neuroanatomical basis of non-human primate 
prefrontal-cerebellar networks is relatively well established, and 
becoming clearer in rodents (see above), scant information is 
available on the dynamic physiological interactions between the 
two structures, particularly during behavior. This is of impor- 
tance given that temporally organized, large-scale, distributed 
networks are thought to be fundamental to information process- 
ing as reflected in frequency specific, coherent local field potential 
(LFP) oscillations (Gray, 1994; Varela et al, 2001; Fries, 2005). 
Indeed, coherent cerebro-cerebellar oscillations have been pre- 
viously observed, across frequencies ranging from 1 to 40 Hz, 
in both anesthetized and awake, behaving animals (O'Connor 
et al, 2002; Courtemanche and Lamarre, 2005; Soteropoulos and 
Baker, 2006; Ros et al, 2009; Rowland et al., 2010). Nevertheless, 
it remains unknown whether prefrontal cortical activity also 
synchronizes with the cerebellum. 

We therefore sought to address the basis and nature of 
cerebello-prefrontal interactions in rat by: (a) using electro- 
physiological mapping techniques to study connectivity between 
the cerebellar vermal output nucleus, fastigius and neocorti- 
cal regions, including the PrL; (b) exploring the possibility that 
cerebellar stimulation may modulate ongoing firing patterns in 
neocortical regions; and (c) examining the coordination of LFP 
activity within this network. 

METHODS 

All experimental procedures were carried out in accordance 
with the UK Animals (Scientific Procedures) Act 1986 and were 
approved by the University of Bristol institutional animal licence 
advisory group. A total of 13 adult rats were used in two 
experimental groups: non-recovery electrophysiology (8 Wistar 
rats, weight 280-380 g, Harlan, UK) and chronic, recovery elec- 
trophysiology (4 Long-Evans and 1 Wistar, weight 340-440 g, 
Harlan, UK). Strain-related differences were not apparent in any 
of the results. 

ELECTROPHYSIOLOGY IN ANESTHETIZED RATS 

Rats were anesthetized with urethane (1.5g/kg intraperitoneal 
injection) then placed in a stereotaxic frame (David Kopf instru- 
ments, Tujunga, CA) and secured with atraumatic ear bars coated 
with a topical local anesthetic (Xylocaine, Astra Pharmaceuticals, 
Kings Langley, UK). Occasionally a supplementary dose of ure- 
thane was given (10 % of original dose) to maintain surgical levels 
of anesthesia, as evidenced by the absence of limb withdrawal 
and corneal reflexes and lack of whisking. Core body temperature 
was maintained at 36-38°C through the use of a homoeothermic 
blanket (Harvard apparatus, Massachusetts, USA). Craniotomies 
were made over the frontal cortex (+3.2 mm, +0.6 mm from 
bregma) and cerebellum (—11.5mm, +0.8 mm from bregma). 

Cortical recordings were made using a Cheetah 32 system 
(Neuralynx, Montana, USA), with extracellular action potentials 
(sampled at 32 kHz and filtered between 0.6-6 kHz) recorded dif- 
ferentially using a local reference placed in a proximal cortical 
region in which spiking activity was absent. Typically, arrays of six 
extracellular tetrode recording electrodes were positioned on the 
surface of secondary motor cortex (M2) whilst a bipolar stimulat- 
ing electrode (SNE-100X, interpolar distance of 0.5 mm, Rhodes 
Electromedical) was targeted toward the contralateral cerebellar 



fastigial nucleus (FN, 4.5 mm from surface of brain) and used 
to deliver trains of stimuli at 0.03 Hz (intensity of 100 u,A and 
frequency of 100 Hz; cf. Mittleman et al., 2008). Extracellular 
responses to the cerebellar stimulation were recorded at two 
depths within the frontal cortex: firstly in superficial regions 
(M2/anterior cingulate border; 1.3-2 mm ventral from brain sur- 
face) and once again when the tetrodes had been lowered to their 
final position in the PrL region (2.6-3 mm ventral from brain sur- 
face). On average, we recorded 2.4 ± 0.2 well-isolated units per 
tetrode in anesthetized animals; continuous LFP was not recorded 
in these anesthetized preparations. 

ELECTROPHYSIOLOGY IN CHRONICALLY IMPLANTED RATS 

Rats were implanted with up to 8 tetrode recording electrodes into 
the left frontal cortex (+3.2 mm, +0.6 mm from bregma) and 
1 bipolar stimulating/recording electrode into the contralateral 
FN (—11.5 mm, +0.8 mm from bregma; interpolar distance of 
0.5 mm) under sodium pentobarbital recovery anesthesia. In one 
animal, tetrodes were implanted in both cerebellum and frontal 
cortex at the same coordinates as given above. Following surgery, 
the independently moveable tetrodes were lowered into the PrL 
subdivision of the prefrontal cortex (~2.6— 3 mm ventral from 
brain surface) over the course of 1 week. Differential recordings 
were made using a Digital Lynx system (Neuralynx, Montana, 
USA) with a local reference placed in a proximal cortical region 
without spiking activity (2.4-2.7 mm below the pial surface for 
prefrontal recordings). On average we recorded 2.2 ± 0.6 units 
per tetrode in the chronically implanted rats. 

Cerebellar LFP recordings were made through either previ- 
ously implanted bipolar electrodes positioned in the cerebellum, 
with overlying skull screws serving as the reference point or in one 
case with tetrodes, which were referenced locally to a tetrode with- 
out spiking activity. LFP signals were sampled at 2 KHz and fil- 
tered between 0. 1 and 475 Hz. Extracellular action potentials were 
sampled and filtered as for the acute, non-recovery experiments 
and recording channels were grounded to two screws in the skull 
overlying the cerebellum. In some cases electrodes were coated in 
Dil (l,l'-dioctadecyl-3,3,3',3'-tetramethylindocarbocyanine per- 
chlorate; Molecular Probes, Invitrogen, UK) prior to implanta- 
tion, which, in addition to electrolytic lesioning, was used to help 
establish electrode tip positions at the end of each experiment (see 
Figure 1). 

CHRONIC RECORDING AND STIMULATION PROTOCOLS 

Evoked field potentials and single unit responses were recorded in 
the M2/anterior cingulate (Cgl) and PrL regions whilst animals 
were in a rest box, which consisted of an elevated platform (20 cm 
diameter) inside a wooden box (45 x 45 x 100 cm). In all rest 
box experiments animal movement was monitored continuously 
by video. For field potential experiments, stimulation parameters 
consisted of a triplet burst of 3 pulses (0. 1 ms pulse duration, 3 ms 
inter-pulse interval) delivered once every 2 s. The mean stimula- 
tion intensity required to evoke reliably detectable field potentials 
was 300 ± 115 [iA (range 100 uA to 500 uA; n = 4). For single 
unit experiments, trains of stimuli (100 Hz, 100 stimuli, 1 s dura- 
tion, mean stimulus intensity 80 ± 20 (xA, range 40 [iA to 100 [iA; 
n = 3) were delivered to FN every 5 s. These stimulus parame- 
ters have previously been shown to drive cerebellar nuclear output 
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FIGURE 1 | Extracellular tetrode recordings of single-unit PrL activity 
during cerebellar stimulation. (A) Grouped schematic (shaded gray areas, 
adapted from Paxinos and Watson, 2006) and representative micrographs 
of neutral red-stained 50 (im saggital and transverse brain slices showing 
respectively the sites of electrolytic lesions (arrowheads) in the cerebellum 
(left,) and PrL, right). Dashed lines indicate outline of fastigial nucleus and 
overlying cerebellar cortical lobules are numbered; scale bars, 1 mm. M2, 
supplementary motor cortex; Cgl, cingulate cortex; IL, infralimbic cortex. 



FMI, forceps minor of the corpus callosum. (B) Three clusters of action 
potentials spread along the axes of relative energy recorded on two 
channels of a tetrode in PrL. The properties of the black cluster (circled) 
are shown in (CD). (C) Mean waveform recorded on color-coded channels 
of the tetrode (top) showing stable relative spike amplitudes throughout 
one experiment (bottom) (scale bar, 0.05 mV; 0.7 ms). (D), Distribution of 
interspike intervals (ISI) for all spikes fired by the unit in the experimental 
session. 



(Bagnall et al., 2009). For experiments in which recordings were 
made from both cerebellum and frontal cortex (n = 5), rat loca- 
tion was video tracked in the open field (aim diameter circular 
arena) via light- emitting diodes attached to a powered headstage 
(Cheetah software; Neuralynx, Montana, USA). 

DATA ANALYSIS 

All data were processed in Matlab (Mathworks, USA) unless 
stated otherwise. LFP and single unit data were sorted based 
upon running speed (derived from video tracking data). For open 
field, a thresholding algorithm extracted stretches of LFP that fell 
within periods of active locomotion, defined as a z-score nor- 
malized running speed of greater than 0. These LFP sections 
were then used for subsequent analysis. For rest box record- 
ings, LFP was selected when the rats were in a state of quiet 
wakefulness characterized by minimal locomotion and absence of 



frontal cortical sleep-spindle activity. Multitaper spectral analyses 
were performed using the Chronux toolbox (Bokil et al., 2010). 
Directed coherence — which uses autoregressive models of two 
LFP signals to estimate which signal best predicts the other — was 
calculated using custom scripts described and published else- 
where (e.g., Baker et al., 2006; Williams et al., 2009). Single units 
were manually isolated off-line (Figures 1B,C) using clustering 
software (MClust3.5; A.D.Redish, available at http://redishlab. 
neuroscience.umn.edu/MClust/MClust.html); inclusion criteria 
were set to isolation distance > 15.0 and L-ratio <0.35 (cf. Harris 
et al., 2001). Putative pyramidal cells were classified on the basis 
of spike width, waveform and mean firing rate (Jung et al., 1998). 
Cross-correlograms were computed in Matlab and spike trains 
shuffle-corrected across trials then normalized by the number of 
spikes. Autocorrelograms were constructed in the same manner 
and normalized by the number of spikes. 
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Peri-stimulus histogram (PSTH) plots were calculated for 
2 s pre- and 5 s post-stimulation epochs with mean baseline 
firing rate calculated from the pre-stimulus period. Firing rates 
were computed in 100 ms bins ± bootstrapped error estimate. 
Trial-averaged rate was calculated and smoothed by a Gaussian 
kernel. The standard deviation of the kernel was set to 0. 1 s. 
Significance in firing rates was determined with a random per- 
mutation test performed with a minimum of 10,000 random- 
izations (cf. Hagan et al., 2012). Significance was assumed when 
mean ± bootstrapped error estimate was above/below the pre- 
stimulation baseline firing rate. Cell response characteristics were 
calculated using automated Matlab scripts and compared using 
X 2 -tests and One-Way AN OVA with post-hoc Tukey's Multiple 
Comparison Test. 

Evoked field potential data were taken from a single tetrode 
channel and averages created in Spike2 software (Cambridge 
Electronic Design, UK). Evoked field potential onset latencies 
were measured from the time of the third and final stimulus of 
stimulation bursts to avoid contamination by stimulus artefacts 
(since total burst duration was 6 ms; see Figure 2A). Since abso- 
lute field potential amplitudes varied between animals, results 
were normalized and expressed as a percentage of the maxi- 
mal response size. Response averages were then compared using 
One-Way ANOVA with post-hoc Tukey's Multiple Comparison 
Test. Data are presented as mean ± s.e.m. unless stated other- 
wise. 

RESULTS 

CEREBELL0-PREFR0NTAL CONNECTIVITY DURING QUIET REST 

We first examined the effects of FN stimulation on LFP in the 
frontal cortex in rats during quiet rest. In three out of four experi- 
ments FN stimulation evoked the largest field potentials at depths 
of 2.6-3 mm from the cortical surface. This depth range corre- 
sponds to PrL (see Figure 2). The evoked field potentials had an 
average onset latency of 13.1 ± 1.1 ms and peak-to-trough ampli- 
tude of 0.23 ± 0.12 mV. By comparison, superficially positioned 
tetrodes recorded evoked field potentials that were ~60% smaller 



FIGURE 2 | Evoked field potentials in the frontal cortex following FN 
stimulation in behaving rats. (A) Example experiment illustrating averaged 
(thick black line) field potentials (72 trials) recorded from tetrodes at different 
depths in the frontal cortex following stimulation of the FN (recording 
positions indicated by numbers on rat brain schematic adapted from Paxinos 
and Watson (2006); small arrowheads indicate timing of FN stimulation 
artefacts; M2, supplementary motor cortex; Cg1, cingulate cortex; IL, 



in size (0.09 ± 0.03 mV; onset latency of 13.2 ± 1.5 ms Figure 2B). 
Recordings from ventral PrL (depth 3.1-3.5 mm) revealed evoked 
field potentials that were ~30% smaller than those recorded at 
2.6-3 mm (0.16 ± 0.04 mV; onset latency of 13.2 ± 1.6 ms). 
Although changes in size of LFPs in the cerebral cortex should 
be interpreted with caution because current source-sink relation- 
ships are complex, nonetheless, the systematic variation in field 
potential amplitude found in the present study raises the pos- 
sibility that this reflects a preferential physiological connectivity 
between FN and PrL compared to other areas of frontal cortex 
that were sampled. 

MODULATION OF PREFRONTAL FIRING FOLLOWING FN STIMULATION 

Following the discovery of cerebellar-prefrontal connectivity at 
the field potential level, we next sought to examine whether 
FN stimulation could modulate the ongoing firing patterns of 
individual frontal cortical neurons. 

The large projection neurons in the cerebellar nuclei can be 
driven to fire at > 100 Hz in slice preparations (Bagnall et al, 
2009). Therefore, as an initial step we examined the effect of 
high frequency FN stimulation (100 Hz, 100 stimuli, Is dura- 
tion, 100 |xA) on PrL cell firing rates in the awake, behaving rat. 
Recordings were made from 20 cells in 3 animals as they sat qui- 
etly on an elevated rest platform (see methods for further details) 
whilst the contralateral FN was stimulated. Of the cells recorded, 
50% (10 cells) displayed a decreased firing rate compared to 
baseline activity (mean baseline firing rate = 5.8 ± 1.2 Hz, see 
Table 1), 5% (1 cell) showed a significant increase, whereas 40% 
(8 cells) displayed a biphasic response (see Figure 3). Only 5% of 
the sample (1 cell) exhibited no change in firing rate following 
stimulation. 

Consistently rhythmic cell firing was not detected in the awake 
animal, and FN stimulation did not modulate either auto- or 
cross-correlations (data not shown), most likely due to stochastic, 
behavior-dependent PrL firing in the awake, behaving rat (Jung 
et al., 1998). Therefore, in order to examine the influence of cere- 
bellar stimulation on more stationary frontal cortical firing, we 




1.3-2 2.6-3 3.1-3.5 

Recording Depth (mm) 

infralimbic cortex. FMI, forceps minor of the corpus callosum); scale bars, 
1 mm and 0.1 mV, 20 ms, respectively. Thin gray lines idicates s.e.m. Arrow 
indicates field potential peak. (B) Grouped field potential peak-to-trough 
amplitudes expressed as a percentage of the maximal reponse size at 
superficial (1.3-2 mm) and intermediate (2.6-3 mm) and ventral (3.1-3.5 mm) 
recording positions (*P < 0.05; One-Way ANOVA with Tukey's multiple 
comparison test; n = 4; each data point calculated from 72 trials per animal). 
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Table 1 | Frontal cortex cell responses following FN stimulation. 



Cell 


# 


Baseline 


% of cells 


Max 




Time of 




% cells 


Minimum 


Time of 


% cells 


% 


recording 


cells 


firing 




increasing 


firing 




max firing 




decreasing 


firing rate 


minimum 


showing 


unresponsive 


depth/ 




rate 




firing only 


rate 




rate (s) 




firing only 


(%) 




firing rate 


biphasic 


cells 


anaesthetic 




(Hz) 






(%) 
















(s) 


response 




<2.5mm 
































(urcthane) 


55 


1.58 ± 
0.20 




28 


460.3 ± 
60.8 




0.66 ± 
0.07s 




17 




29.3 ±- 
1.8 




L.85± 
0.15 


55 


0 


2.6-3mm 
(urethane) 


69 


1.10 ±" 
0.12 




18 


497.4 ± 
47.9 


* 

* 


1.28 ± 
0.15 




42 


* 
# 


1 12.4 ± 
2.3 


• * 

* 
* 


1.62 ±. 
0.13 


33 


7 


2.6-3mm 
(awake) 


20 


5.8 ± 
1.2 




5 


148.8 
±12.0 




1.63 ± 
0.33 


50 




50 ± 
6.4 


0.94 ± 
0.2 J 


40 


5 



P-values between data pairs marked by brackets. One-Way ANOVA with Tukey's multiple comparison test between cells recorded at <2.5mm depth and those at 
2.6-3 mm under anesthesia, and cells recorded between 2.6 and 3 mm in anesthesia and awake states. *P < 0.05; **P < 0.07; ***p < 0.007. 



next examined the result of FN stimulation on putative pyra- 
midal cell firing rates in rats anesthetized with urethane, which 
induces dominant slow-wave oscillations (cf. Clement et ai, 
2008) that have been suggested to play an important role in 
neocortical-cerebellar communication (Ros et al, 2009; Rowland 
etal., 2010). 

Using the same parameters as for awake rats (100 Hz, 100 
stimuli, 1 s duration, 100 |iA), we found that stimulation of 
the contralateral FN resulted in a robust but heterogeneous 
modulation of frontal cortex putative pyramidal cell firing (see 
Figure 3 and Table 1 for comparison of response characteris- 
tics). From a total of 55 cells recorded (n = 5 rats) in M2 and 
Cgl regions of frontal cortex under urethane anesthesia, 27.3% 
(15 cells) displayed a significant increase in firing rate com- 
pared to pre-stimulation baseline activity, whereas 16.4% (9 cells) 
showed a significant firing rate decrease. A third category of 
cells (56.3%; 31 cells) showed a biphasic response that gener- 
ally consisted of a significant firing rate decrease and subsequent 
increase (see Table 1). Cells recorded from tetrodes positioned 
at depths corresponding to PrL [n = 8 rats; 69 cells) were also 
heterogeneously modulated by the FN stimulation but the over- 
all pattern of responses observed in PrL cells was significantly 
different to those recorded in M2/Cgl (PrL cells: 12/69 (17.4%) 
displayed an increase in firing rate, 29/69 (42 %) displayed a 
decrease in firing rate, 23/69 (33.3%) showed a biphasic pat- 
tern, and 5/69 (7.3%) showed no response (Figure 3C, when the 
proportion of cells in these different categories were compared 
to those found in M2/Cgl x 2 = 15.66, df = 3, **P < 0.01; see 
Table 1). 

Of the cells that responded to FN stimulation with decreases in 
firing rates, PrL cells showed more profound firing rate reductions 
than M2/Cgl cells (to 29.3 ± 1.8% and 12.4 ± 2.3% of base- 
line respectively; ***P < 0.001, One-Way ANOVA with Tukey's 



multiple comparison test; see Table 1). Of cells responding with 
increased firing rates, relative increases were similar in PrL 
and M2/Cgl populations (to 497 ± 47.9% and 460 ± 60.8% 
respectively), though FN stimulation-induced firing peaked more 
rapidly in M2/Cgl than in PrL (0.66 ± 0.07 s following stimu- 
lation compared to 1.28 ± 0.15 s **P < 0.01, One-Way ANOVA 
with Tukey's multiple comparison test). Overall, these data there- 
fore reflect a complex pattern of modulation, with M2/Cgl cells 
tending to respond with a rapid biphasic response and PrL cells 
typically displaying an initial reduction in firing rate following 
cerebellar stimulation. 

This overall pattern of response did not differ significantly 
from the equivalent recordings made in awake rats (x 2 = 2.2, 
df = 3, P > 0.05; see Table 1). However, compared to recordings 
made in awake rats, PrL cell firing rate increases/decreases in anes- 
thetized animals were significantly more pronounced following 
FN stimulation (See Figure 3 and Table 1). 

Next, by using auto- and cross- correlogram analyses, we inves- 
tigated the effect of FN stimulation on the average, coordinated 
network rhythmicity within the PrL in urethane anesthetized 
rats. Despite the heterogeneity in PrL cell responses shown in 
Figure 3 and Table 1, cerebellar stimulation resulted in modu- 
lation of ongoing population PrL network activity, as observed 
in the disruption of slow wave oscillations, and broadening 
of the central peak in both auto- and cross-correlations (see 
Figure 4). This finding highlights the potential of the cerebellum 
to influence ongoing network processing in the neocortex and 
provides further evidence of functional connectivity between the 
regions. 

CEREBELL0- PREFRONTAL COMMUNICATION IN AWAKE RATS 

As a first step to understanding cerebello-prefrontal network 
activity and interactions in behaving animals, we examined 
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FIGURE 3 | Single unit PrL responses following FN stimulation in 
awake and urethane anesthetized rats. Raster and peri-stimulus rate 
plots for example cells recorded in PrL in the awake animal, cells 1 
and 2 (A) and urethane anesthetized animal, cells 3 and 4 (B). 
Horizontal black bar indicates duration of stimulation (100(aA; 100Hz; Is 
duration); bold line indicates instantaneous mean firing rate; gray lines 



indicate bootstrapped error estimate; horizontal line indicates mean 
baseline firing rate prestimulation; bin size, 100 ms; 17 trials for 
anesthetized and 33 trials for awake experiments. (C) Quantification of 
PrL cell firing patterns following FN stimulation in awake (n = 20 
putative pyramidal cells from 3 animals) and anesthetized rats (n = 69 
putative pyramidal cells from 8 animals). 



the covariance of prelimbic cortical and fastigial nucleus 
LFP signals using Fourier coherence analysis (see meth- 
ods for further details) during active locomotion in a 1 m 
diameter open field arena vs. quiet restfulness on a 20 cm 
platform. 

Cerebellar theta power (5-10 Hz) showed a slight increase dur- 
ing active locomotion relative to rest (Figure 5A top; FN theta 
power during rest, 17 ± 1.6 dB; active locomotion, 25 ± 2.3 dB; 
P < 0.05, Wilcoxon rank sum test; n = 4 for open field and 5 
for rest box recordings), whilst PrL theta power was similar in 
the two behavioral states (PrL theta during rest, 30 ± 4.8 dB; 
active locomotion, 29 ± 6.0 dB; P > 0.05). Despite these lim- 
ited power changes, the FN LFP signal was significantly and 
selectively coherent with PrL oscillations in the theta range (5- 
10 Hz) only during active locomotion in the open field (P < 0.05, 
arrow in Figure 5B), but not while rats were at rest (Figure 5A). 
Consequently, the proportion of total coherence carried at theta 



frequency (a ratio between 5-10 Hz coherence and coherence at 
all other frequencies up to 45 Hz) was significantly higher during 
locomotion (ratio 1.85 ± 0.21) than rest (0.97 ± 0.01; P < 0.05, 
Wilcoxon rank sum). 

Since coherence is a measure of consistent phase relation- 
ships and does not quantify the direction of interaction between 
two signals, we also used Granger causality (directed coher- 
ence; see methods for details) to infer directionality from the 
simultaneous FN and PrL LFP recordings (Figure 5). Uni- 
directional theta coherence was significantly weighted in the 
FN-PrL direction when animals were actively moving in the 
open field (see Figure 5B lower panel; P < 0.05 FN-PrL vs. 
PrL-FN, Wilcoxon rank sum test; n = 4). In contrast, theta 
coherence was no longer significantly directional (FN-PrL and 
PrL-FN comparision P > 0.05) and was significantly lower dur- 
ing rest (P < 0.05 vs. active locomotion, Wilcoxon rank sum 
test). 
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These coherence analyses provide evidence to suggest that the 
cerebellar-prefrontal connectivity exemplified by FN stimulation- 
evoked responses in PFC could subserve cross-structural, network 
interactions that preferentially manifest at theta frequencies dur- 
ing active behavior. 



DISCUSSION 

The key findings of the current study are that: (1) stimulation 
of the FN evokes short latency (~13ms) field potentials in PrL 
alongside changes in frontal cortical neuronal firing rates and 
rhythmicity and (2) FN and PrL LFP are coherent in the theta 
(5-10 Hz) frequency range during active locomotion, an effect 
preferentially driven in the FN-PrL direction. These findings 
demonstrate physiological interactions between vermal cerebel- 
lum and prelimbic cortex in rat and provide insights into the 
neural dynamics of the reciprocally connected networks under- 
pinning cerebellar-cerebro communication. 

CEREBELLAR-PREFRONTAL CONNECTIVITY 

Recent physiological evidence from rats indicates the presence of 
a prefrontal-olivo-cerebellar projection specifically to vermal lob- 
ule VII (Watson et al, 2009). In the current study we targeted the 
fastigial nucleus, the output of vermal lobule VII known to inte- 
grate signals from the cerebellar cortical A-zone (e.g., Voogd and 
Ruigrok, 2004), to examine reciprocal cerebellofugal influence 
on the prefrontal cortex. Our finding that fastigial stimulation 
can evoke discrete field potentials in the PrL (Figure 2) cor- 
roborates previous electrophysiological studies: fastigial nucleus 
projects to widespread cerebral cortical areas via the ventrome- 
dial thalamic nuclei in cat (~10 ms latency; Steriade, 1995), plus 
cerebellar dentate nucleus stimulation results in short latency field 
potential responses in the prefrontal association areas of monkeys 
(2-4.5 ms latency; Sasaki et al, 1979). 

Although we cannot categorically rule out incidental stimu- 
lation of neighboring cerebellar nuclei, our histological verifi- 
cation of stimulation sites, inter-animal consistency of LFP and 
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FIGURE 4 | Slow-wave modulation of PrL activity is disrupted by FN 
stimulation in urethane anesthetized rats. Cross- and auto-correlogram 
plots (40 ms bins, calculated over 2 s pre/post-stimulation) of all possible 
PrL cell pair combinations (n = 69 cells, 8 rats) during non-stimulated (thin 
line) and FN stimulation (thick black line) states in urethane anesthetized 
rats. **P < 0.01, paired f-test. Note attenuation of slow-wave periodicity 
during FN stimulation also reported in anesthetized cat (Steriade, 1995). 




FIGURE 5 | PrL-FN LFP network activity during rest and open field 
exploration. Grouped power spectra (top panel; cerebellum/FN in blue, 
prelimbic cortex/PrL in black), coherence (middle panel) and directed 
coherence (Dir coh, bottom panel; FN-PrL direction in blue, PrL-FN in 
black) as rats sat quietly on rest platform (A, n=5) or actively moved 



in the open field arena (B, 1 Hz bandwidth; n = 4). Confidence level at 
P = 0.05 marked by horizontal black lines, shading indicates jack-knife 
error bars. Arrow indicates significant coherence peak at theta 
frequency, which was only during active locomotion and preferentially 
driven in the FN-PrL direction. 
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comparable results in cat lend weight to the fastigial nucleus 
constituting a key relay in cerebello-cerebral interactions. The 
latency observed in the current study (~13 ms) suggests that the 
cerebello-prefrontal pathway in rat is relatively slow conducting 
and/or polysynaptic, presumably involving at least one synaptic 
relay (most likely within the thalamus), though it is possible that 
a direct pathway to a frontal subregion other than M2/Cgl/PrL 
does exist. Alternatively, the shorter latency and faster conduction 
of the cerebello-prefrontal pathway in non-human primate could 
reflect evolution of a more rapid, direct line of communication in 
line with the increasing importance of this pathway in cerebellar 
contributions to cognitive functions (Ramnani, 2006). 

CEREBELLAR MODULATION OF PREFRONTAL FIRING 

Electrical stimulation of the fastigial nuclei drives heterogeneous 
modulation of prefrontal cell firing in both anesthetized and 
behaving rats. Preliminary findings in awake mice have reported 
similar frontal responses to cerebellar cortical stimulation, show- 
ing reductions in synchronized firing and pauses in cortical cell 
firing (Liu et al, 2012). Steriade (1995) showed that fastigial 
stimulation can modulate gamma frequency (~20-40Hz) EEG 
rhythms in the frontal cortex of anesthetized cats (with effects 
outlasting the duration of stimulus by ~4-5 s) and cerebellar 
cortical stimulation has recently been shown to modulate EEG 
activity recorded over the motor cortex of awake mice (Ml, 
Witter et al., 2013), thus highlighting the presence of multiple 
cerebello-cortical pathways. 

Mittleman et al. (2008) described long lasting (7-8 s) 
dopamine efflux in the PrL of mice following cerebellar den- 
tate nuclear stimulation and dopamine release may contribute 
to the effects on PrL found in the present experiments. Primate 
dorsolateral prefrontal cortex pyramidal neurons express mixed 
dopamine receptor distributions (Lidow et al., 1991; Muly et al., 
1998), and in vitro experiments have highlighted the opposing 
effects that Dl and D2 receptor types exert on prefrontal neuron 
spiking (Seamans et al., 2001). The heterogeneous effects of fasti- 
gial nuclear stimulation on individual pyramidal cell responses 
summarized in Table 1 may therefore reflect mixed dopamine 
receptor subtype expression and/or targeting by dopaminergic 
afferents. 

The pathway(s) through which cerebellar stimulation can 
influence prefrontal dopamine release are unknown but one 
potential route includes projections via the thalamic nuclei, which 
send glutamatergic afferents to the cerebral cortex (Hoover and 
Vertes, 2007) that in turn form presynaptic inputs to dopamine 
varicosities within the prefrontal cortex (leading to a slow, neu- 
romodulatory response, Blaha et al., 1997). Alternatively, fastigial 
stimulation could influence the release of dopamine in the pre- 
frontal cortex via activation of the ventral tegmental area (VTA) 
through cerebello-VTA or cerebello-thalamo-VTA projections 
(Kehr et al, 1976; Snider and Maiti, 1976; Snider et al, 1976). 
The cerebellum may also influence prefrontal firing via the basal 
ganglia. For example, electrical stimulation of the cerebellar out- 
put nuclei (dentate in particular) alters neuronal firing rates (Li 
and Parker, 1969; Ratcheson and Li, 1969) and dopamine levels in 
the substantia nigra and caudate nucleus (Nieoullon et al, 1978). 
Anatomical studies demonstrate disynaptic projections from the 



cerebellum (dentate) to the basal ganglia in both rat and monkey 
(Ichinohe et al., 2000; Hoshi et al., 2005). In turn, basal gan- 
glia projections influence prefrontal cortex (e.g., Middleton and 
Strick, 1994; Maurice et al, 1999; Middleton and Strick, 2002), 
thus highlighting the basal ganglia as a potential relay between 
hind- and forebrain. 

COHERENT PREFRONTAL-CEREBELLAR LFP DURING BEHAVIOR 

Neurobiological oscillations organize activity within and across 
different brain regions (e.g., Singer, 1999; Varela et al., 2001; Fries, 
2005), creating coherent cell assemblies (Harris et al., 2003) and 
enabling plasticity processes dependent on the precise timing of 
pre- and post-synaptic activity (Markram et al., 1997; Bi and Poo, 
1998; Cassenaer and Laurent, 2007). Coordinated oscillations 
may therefore support information transfer in cerebro-cerebellar 
pathways and have been reported across a range of frequencies: 
cerebellar oscillations phase-lock to neocortical slow waves (0- 
4 Hz; Ros et al, 2009; Rowland et al., 2010; Schwarz, 2010) and 
beta oscillations (~ 10-25 Hz) bind cerebellar cortical LFP and 
nuclear cell firing with the somatosensory and motor cortices of 
primates (Courtemanche and Lamarre, 2005; Soteropoulos and 
Baker, 2006). 

Of particular relevance to the current study are (1) the 
demonstration by Steriade (1995) that stimulation of cat FN dis- 
rupts slow wave activity, a result corroborated by our findings 
(Figure 2) and (2) the additional finding that cerebellar circuits 
have been shown to support oscillation frequencies within the 
theta bandwidth (Hartmann and Bower, 1998; D'Angelo et al., 
2001; Solinas et al, 2007; Dugue et al, 2009), which could 
be driven by pacemaker theta rhythmicity within precerebellar 
nuclei including the inferior olive (Lang et al, 2006; Chorev et al., 
2007; Van Der Giessen et al., 2008). Theta rhythmicity in the cere- 
bellar fastigial nuclei may therefore reflect synchronous activation 
or inhibition in either the olivo-nuclear or olivo-cerebello-nuclear 
circuitry, and there is evidence that such oscillations synchronize 
activity within cerebellar hemispheres as well as between cere- 
bellar and cortical/limbic regions (Hartmann and Bower, 1998; 
O'Connor et al., 2002; Hoffmann and Berry, 2009; Wikgren et al., 
2010). 

O'Connor et al. (2002) found that LFP activity in whisker- 
related areas of rat cerebellar cortex and neocortex are coherent 
at 5-20 Hz during periods of active whisking, and whisker- and 
eye-movement related activity is found in the rat prefrontal 
region (Brecht et al., 2004). Coherent cerebellar-cerebro activ- 
ity may therefore reflect mechanisms through which sensory 
information can be integrated into ongoing neocortical pro- 
cesses (cf. Bland and Oddie, 2001; O'Connor et al, 2002; Bland, 
2004). Our finding that coherence is significantly weighted in the 
cerebellum-to-PrL direction during epochs of active locomotion 
(Figure 5B) suggests the interaction derives from more than sim- 
ple co-modulation of cerebellum and PrL during theta-frequency 
behaviors (e.g., whisking). The directed nature of this coupling 
may reflect the increased need for sensorimotor input to the 
neocortex during goal-directed behaviors including active loco- 
motion. In particular, as vestibular information is combined with 
proprioceptive inputs in the FN to generate appropriate inter- 
nal estimates of the animal's self motion (Brooks and Cullen, 
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2009), the cerebellum (FN in particular) may provide function- 
ally relevant proprioceptive/egocentric information that can be 
integrated into decision making processes recruiting higher order 
structures such as the PrL. 

CONCLUDING COMMENTS 

Cerebellar vermal abnormalities are found in a host of psychiatric 
diseases (Heath et al, 1979, 1982; Okugawa et al, 2007; Lawyer 
et al., 2009) and can occur concomitantly with changes in the 
prefrontal cortex of autistic patients (Carper and Courchesne, 
2000). Also, chronic cerebellar vermal stimulation in the theta 
frequency range has been reported to ameliorate the emotional 
and cognitive symptoms of intractable neurological disorders 
such as schizophrenia and epilepsy (Cooper, 1973; Cooper et al., 
1976; Correa et al, 1980). Abnormalities of the cerebellum, 
and particularly its vermal region, may therefore contribute to 
neuro-psychiatric diseases that are typically associated with neo- 
cortical malfunction and aberrant dopamine neuromodulation. 
The present results provide evidence for a physiological frame- 
work whose dysfunction could underlie cerebellar contributions 
to such disorders. 

Further work monitoring and manipulating cerebello-cerebral 
network activity with higher resolution during a range of behav- 
ioral states (e.g., with and without explicit cognitive load) is 
required before we fully appreciate the functional importance 
of activity in cerebello-prefrontal circuits. However, the current 
study provides initial evidence that the regions co-participate 
in distributed network activity and also offers novel insights 
into the dynamic interaction that occurs between the two struc- 
tures during exploratory behavior. The potential mechanisms 
subserving these interactions include synchronized oscillations 
in the theta frequency, which may be important for sensory 
acquisition/integration (Bower, 1997; Bland and Oddie, 2001) 
and/or general cerebellar contributions to goal directed behaviors 
(Burguiere et al., 2010). 
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